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ABSTRACT 

Proteins from the kinesin-8 family promote microtubule (MT) depolymerization, a process thought 
to be important for the control of microtubule length in living cells. In addition to this MT short- 
ening activity, kinesin 8s are motors that show plus-end directed motility on MTs. Here we describe 
a simple model that incorporates directional motion and destabilization of the MT plus end by ki- 
nesin 8. Our model quantitatively reproduces the key features of length-vs-time traces for stabilized 
MTs in the presence of purified kinesin 8, including length-dependent depolymerization. Compari- 
son of model predictions with experiments suggests that kinesin 8 depolymerizes processively, i.e., 
one motor can remove multiple tubulin dimers from a stabilized MT. Fluctuations in MT length 
as a function of time are related to depolymerization processivity. We have also determined the 
parameter regime in which the rate of MT depolymerization is length dependent: length-dependent 
depolymerization occurs only when MTs are sufficiently short; this crossover is sensitive to the bulk 
motor concentration. 
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Introduction 

Regulation of microtubule (MT) length is an important cellular process. Abnormal MT lengths 
can mislocalize the nucleus or mitotic spindle and cause defects in polarized cell growth or mitosis. 
Although MT length regulation is not fully understood, several mechanisms have been proposed, 
including chemical gradients in the mitotic spindle jl|), cortical interactions trafficking of pro- 
teins that bind to MT ends (0), and a balance between proteins that promote MT polymerization 
and depolymerization (0, IE) • 

Kinesin-8 proteins appear to help regulate MT length in vivo. Deletion of kinesin-8 genes leads 
to longer interphase and spindle MTs and defects in mitosis (0, 0, 0, 0, Eo , U, H), suggesting that 



kinesin 8s promote MT depolymerization. In addition, kinesin-8 motors show processive, plus-end 



directed motility on MTs in vivo (jl2j, Il3l ). Recent work has shown that kinesin-8 proteins are 
important in chromosome oscillations and MT length fluctuations [l3l. 14). 

Biochemical experiments with some purified kinesin 8s have reproduced many of observations 
made in vivo: kinesin 8 moves processively toward the MT plus end (fla* ) and the MTs then 
depolymerize, even when stabilized with taxol or gmpcpp (filial I la ) . Kinesin 8s are thus thought 



to be able to remove the GTP-tubulin cap that stabilizes growing MTs in vivo. Varga et al. recently 
proposed that the processive motility of kinesin 8s preferentially increases their concentration at 
the ends of longer MTs (fril ). thereby more rapidly depolymerizing longer MTs. They proposed that 
this length-dependent activity serves to regulate MT length. Howard and Hyman proposed that a 
constant MT growth rate coupled with length-dependent depolymerization would result in a tighter 
distribution of MT lengths than that set by dynamic instabilit y M ). In particular, MT dynamic 
instability gives an exponential distribution of MT lengths (0, fl7l)7 while a coupling between MT 
growth and length-dependent depolymerization could give a tighter distribution (fl8l ). The in vitro 
experiments are consistent with the in vivo observation of longer MTs when the kinesin 8 is deleted, 
but the reasons why kinesin 8 depletion and overexpression alter mitotic oscillations are less clear. 

Given that not all of the in vivo results can be simply understood in terms of the in vitro 
observations, we sought to determine if a detailed theory of the in vitro experiments could lead 
to insights into the behavior of kinesin 8s in cells. We have developed a simple mathematical 
model that quantitatively captures the length-vs-time traces of stabilized MTs in the presence of 
purified kinesin 8. The results are consistent with processive MT depolymerization by kinesin 8, i.e., 
multiple tubulin dimers can be sequentially removed by a single motor. We explored a consequence 
of this processive depolymerization: altered fluctuations in MT length during depolymerization. 
In addition, we studied the distribution of motors along the MT and find that length-dependent 
MT depolymerization occurs only for sufficiently short MTs, i.e., those below a "crossover" length, 
which depends strongly on the bulk motor concentration and model parameters. 

Our work builds on previous physical theory that considered the motion of multiple motors on a 
one dimensional track: Parmeggiani et al. studied a model similar to ours but disallowed changes in 
the length of the track (fl^ ) . Nowak et al. extended this work to allow track lengthening catalyzed 
by motors (0). Other theoretical papers have focused specifically on MT depolymerization by 
kinesins. For example, the coupling between motor motion and MT depolymerization is thought to 
be important for kinesin-13 proteins, which also depolymerize GMPCPP- and taxol-stabilized MTs 



(|2ll ). An important difference from the present work is that the kinesin-13 protein MCAK moves on 
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MTs by diffusion along the MT lattice, not through motor activity, and it accumulates at the MT 
ends through a binding preference for this part of the MT. A previously developed theory of MT 
shortening by MCAK focused on the dynamic accumulation of motors at the MT end (|22l ). Both 
the experiments and theory on MCAK found evidence for processive depolymerization. Finally, a 
model of MT shortening catalyzed by either kinesin-8 or -13 motors was recently introduced (|18l ). 

Our work extends previous theories of motor motion and depolymerization in important ways. 
The mathematical model we consider is similar to that of Govindan et al. (18). The key differences 
are that Govindan et al. neglected motor crowding effects and analyzed the steady state of the 
model, while here we include crowding effects and non-steady-state dynamics. An important the- 
oretical challenge arises because the rate at which motors accumulate at the MT end and the rate 
of MT length change are similar (as observed in experiments, see refs. (0, EEI)). A steady-state 
analysis, while mathematically more tractable, may therefore miss dynamic behavior present in the 
experiments. In this paper we first consider a kinetic Monte-Carlo simulation of the full model. 
Then we develop a mean-field model that explicitly considers non-steady-state effects by describing 
the coupling between the dynamics of the MT end and the motor occupancy at the end. The 
depolymerization rates we obtain in the mean-field model agree well with Monte-Carlo simulations 
of the full model. We compare the results of the model, including non-steady-state effects, to the 
experimental data obtained by Varga et al. on length-dependent depolymerization of MTs by 
Kip3p. 

A mathematical model of kinesin-8 action on MTs quantitatively 
reproduces experiments 

Our theoretical model of kinesin-8 motors includes biased motor motion toward MT plus ends and 
motor-catalyzed depolymerization of MT plus ends. A schematic of the model is shown in figure 
[T]A.. The position x, in units of tubulin dimers, is measured from the MT minus end (x = 0). We 
assume motors track individual protofilaments and step from one tubulin dimer to the next (a 
length of 8 nm) toward the MT plus end at rate v. We assume the motors do not step backward 
(&;_ =0), as kinesin-8 motors have shown highly biased motion (jl5l ). 

Motors can bind to and unbind from the MT. The "on" rate is k -n.c, where c is the bulk 
concentration of motors (typically assumed to be constant). The "off" rate is k g, unless the 
motor is bound at the end of the MT. Motors dissociate from the plus end of the MT at rate 
^off^' ^ e ne sl ec t an y special affinity that the motors may have for the MT minus end, since their 
plus-end directed motility makes their occupancy at the minus end negligible. We also assume 
that motors have a negligible affinity for soluble tubulin dimers. Although the binding affinity of 



kinesin-8 motors for soluble tubulin is currently unknown (|12l). when we allowed motor binding to 
free tubulin with an affinity as large as the affinity for MT dimers, there was little change in the 
results (see below, and other data not shown). Therefore, the precise value of the affinity is not 
important for our results. 

We considered two distinct depolymerization mechanisms: one assumes processive depolymer- 
ization (P) and another assumes purely non-processive depolymerization (NP). In the first case (P) 
motors processively depolymerize protofilaments by removing one dimer and stepping backward at 
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rate k^^. If there is another motor bound behind the depolymerizing motor, the motor at the end 
is assumed to unbind. This assumption is motivated by live-cell imaging of fluorescent Kip3p, which 
showed that a clump of motors accumulates at the MT plus end during MT growth, but the clump 
fluorescence greatly decreases during MT shrinkage ( ll la). Since protofilaments are straight in the 
bulk of the MT and protofilament curvature is assodated with depolymerization 0) , our picture 
of processive depolymerization is consistent with a structural mechanism in which the motor has a 
higher affinity for a straight protofilament than a curved region of the protofilament. Our picture 
of processive depolymerization is also consistent with individual motors inducing catastrophe (in 
dynamic MTs) and thereby removing large numbers of tubulin dimers. In the absence of crowding 
effects (if motors fall off before reaching another motor), then the number of dimers removed per 
motor is a = kf^^/k^^. Since motor crowding reduces the depolymerization processivity, the 
value a will only be observed for sufficiently low motor concentration. Note that this picture of 
depolymerization can apply either to a motor that directly catalyzes depolymerization, or a motor 
that recruits to the MT end another protein that promotes depolymerization. 

In the second depolymerization mechanism (NP) motors are assumed to remove a single dimer 
and fall off with that dimer. In this case, each motor removes at most one dimer (a motor could 
unbind before removing a dimer). 

Our idealized MTs consist of 13 protofilaments arranged cylindrically (we neglect helical ar- 
rangement of protofilaments). In all modeling, we further assume that dimer removal by motors 
is non-cooperative, both within and between protofilaments. If each motor acts independently 
(the case of non-cooperative motors), then the rate of depolymerization is, on average, propor- 
tional to the density of motors at the end of the MT. (Note that we do consider the possibility of 
protofilament interactions in tubulin depolymerization; see below.) 

Typically we consider motor- induced depolymerization of stable MTs, so that the MT has no 
intrinsic dynamics. However, in some cases we use a simple model of MT dynamics: dimers add to 
a MT end at rate act (where q is the bulk concentration of tubulin dimers) and dissociate at rate (5 



polymerization rate (3 to depend on the number 
25l). In this case, the rate at which a terminal 



(fig. [T]) . In one set of simulations we allowed the de 
of lateral interactions between protofilaments (24, 
tubulin dimer unbinds from a protofilament is given by (3q if the dimer has no lateral neighbors, (3\ 
if the dimer has 1 lateral neighbor, and if the dimer has 2 lateral neighbors. 

We developed a kinetic Monte Carlo simulation of the full model, and studied analytic approx- 
imations to the model. 



Depolymerization of stabilized MTs 

In the minimal model of stabilized MTs (figure 1A-D), there are 5 independent parameters. Three 
parameters can be derived directly from the data of Varga et al. on the budding yeast kinesin- 
8 motor Kip3p (15). (See table 1 for a summary of estimated parameter values.) For these 
estimates, we assumed that each motor and protofilament behaves independently: we neglected 
effects of protofilament interactions, motor crowding, and motor depolymerization cooperativity. 
The measured motor velocity was 3.6 fim min -1 which gives v = 7.5 dimers s _1 . The average run 
length of a motor was 12 /xm, which gives I = 1500 dimers. (Note that this run length is likely a 
lower bound, since the experiments ignored motors that reach the MT end when determining the 
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run length.) Since I = v/k g, this implies k g = 5 x 10 3 s . The typical residence time of a 
motor at the end of the microtubule was half a minute, so fe^g d = 3.3 x 10 -2 s -1 . 

The two parameters fc^ nc ^ and k on were not directly measured for Kip3p, but we can estimate 
their values from the experimental data. The maximum depolymerization velocity observed was 
approximately 2 /xm min -1 . This suggests k^^ > 4 dimers s -1 . However, by comparing our 
simulations with the data of Varga et al. we found a better fit using the value k^ d = 13 dimers s -1 
(see below). We estimated the motor binding rate from the kymograph of figure 3d from Varga et 
al., which shows 9 binding events in 6 minutes on a MT 12 /im in length, for an unspecified bulk 
motor concentration (of order 1 nM). This gives a microscopic rate constant of order 10 -6 site -1 
nM -1 s -1 . We determined the value k on = 3 x 10 -6 site -1 nM -1 s -1 by comparing simulations 
to the experimental data. In our model, the parameter with the greatest uncertainty is feon- We 
discuss most results using our best-fit value and discuss below the consequences of varying this 
value. 

We first simulated intrinsically stable MTs. The simulations started with bare MTs (no mo- 
tors bound), and motor binding was begun at t = 0. Using the values of k^^ = 13 s -1 and 
k on = 3 x 1(T 6 site -1 nM -1 s -1 , we found good qualitative agreement with experiments. The de- 
polymerization rate in the simulations increased as the motor occupancy along the MT increased. 
The time required for equilibration of the motor occupancy was approximately the typical time a 
motor stayed bound, l/fe ff (see further discussion of occupancy equilibration below in the section 
on motor occupancy profiles). The depolymerization rate later decreased once the MT shortened 
sufficiently 

Our simulations of stabilized MTs did not exhibit the long-time behavior observed in the Varga 
et al. experiments (fl5l ). In these experiments, kinesin-8 motors were not able to fully depolymerize 
GMPCPP-stabilized MTs: the depolymerization rate slowed and dropped to zero over several 
minutes, leading to a "long-time tail" in the plot of MT length as a function of time. This could 
have occurred for several reasons: (i) the motor activity may have decreased, e.g., due to depletion 
of ATP; (ii) motors may have bound to free tubulin dimers or the coverglass, effects which would 



deplete the concentration of free motors (|12l . Il5l ): (Hi) motors may depolymerize cooperatively, so 



that the decrease in motor concentration at the end of a MT due to shortening had a nonlinear effect 
on the depolymerization velocity; or (iv) the slow polymerization activity of GMPCPP tubulin |2^ ) 
may have prevented full MT depolymerization, because tubulin dimers could add to the end of the 
MT, increasing its length. 

We found that motor binding to free tubulin is not likely to explain the long-time tails. We 
simulated binding of motors to free tubulin dimers with varying binding affinity. Even when the 
affinity is large — as large as the affinity for tubulin in the MT lattice — we found only a small change 
in the depolymerization dynamics (data not shown). 

We then modeled polymerization and depolymerization of the MTs with nonzero values of the 
rate constants a and (3 (fig. [T]). The intrinsic dynamics of the MT plus end could have two effects 
on the motors at the MT end: (i) the motors could remain attached to the tubulin dimer that 
is removed or (ii) the motor could be displaced backward and remain on the MT. We found that 
the former led to significantly decreased MT depolymerization rates (the intrinsic MT dynamics 
lead to the unbinding of many motors), while the latter gave quantitative agreement with the in 
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vitro experiments. In particular, we estimated that the total tubulin concentration in the Varga et 
al. experiments is of order 100 nM. Therefore, we included in the model the measured rates a = 5.4 
dimers [iMT 1 s _1 and (3 = 0.1 dimers s _1 (26) and found quantitative agreement between the simu- 
lation traces and the experiments (see figure[2]). In this case, the free tubulin concentration becomes 
high enough that polymerization is large enough to balance motor-induced depolymerization and 
the MTs approach a constant, time-independent length. Thus, we conclude that the accumulation 
of free tubulin within the flow chambers, and the subsequent slow polymerization of the MTs, can 
account for the long-time tails observed in the Varga et al. experiments (fl5l). 

After showing that including the slow polymerization activity of GMPCPP tubulin in the model 
is sufficient to resolve this qualitative disagreement between experiments and theory, we fit the ex- 
perimental data to determine the unknown parameters. With processive depolymerization (model 
P), the only free parameters were fc5. nd , k on , and /3 ,i. Fi gure [3] shows similar MT length versus 
time curves for different values of these parameters. In model P with independent protofilaments 
(A),i = 0) a lower depolymerization processivity could (within some range) be offset by a higher 
on rate (fig. EJ<\). In addition, decreasing the stability of neighborless protofilaments has similar 
effects to increasing the MT depolymerization processivity (fig. [3p). With nonprocessive depoly- 
merization (model NP), only very high /?o,i give MT depolymerization time courses that match the 
experimental results (see below, where nonprocessive depolymerization is discussed in more detail). 
The best-fit experimental traces were obtained with k^ 1 = 13 s~ 1 and fc n = 3 x 10" 6 site" 1 nM" 1 
s- 1 . 



Evidence for processive depolymerization by kinesin-8 motors 

Processive depolymerization of MTs by motors is more consistent with the parameters measured 
by Varga et al. than non-processive depolymerization. In particular, they found that an individual 
motor is resident at the MT end for approximately 30 seconds (at low motor density), and the 
depolymerization velocity they observed had a maximum of 2 ^m min -1 (at higher motor density), 
suggesting that (if residence times are similar over a range of motor densities) approximately 125 
tubulin dimers could be depolymerized during the binding time of a single motor at the MT end 
(fl5l ). Our best fits to the experimental data gave a higher maximum depolymerization rate for 
individual motors. We found the best fit with k^ 1 ^ = 13 dimers s — 1 . This would imply a maximum 
number of dimers removed per motor a = k^^/k^^ ~ 400 dimers. (Note that our simulations 
assumed no interactions between protofilaments or motors, each of which could alter the apparent 
number of dimers removed per motor.) 

We performed simulations to determine whether non-processive depolymerization is consistent 
with the experiments of Varga et al. In this model (NP), we assumed that each motor could remove 
only one tubulin dimer. We tried but failed to match the experimentally observed traces assuming 
completely stable MTs. Increases in the motor on rate by up to a factor of 1000 still produced 
simulations in which depolymerization occurred much more slowly than seen in experiments (data 
not shown). 

Nonprocessive depolymerization by itself is insufficient to account for the experimental results. 
However, combining model NP with intrinsically unstable MTs is partially consistent with exper- 
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iments. The apparent motor processivity increases if protofilaments are not completely stable. 
We altered the protofilament depolymerization rate j3 to depend on the number of neighbors the 
terminal tubulin dimer has. This mode ling choice mimics the stabilizing lateral bonds thought to 
be present between protofilaments (24. l25l). We note that our simple model does not attempt to 
fully describe MT dynamics, and we considered only the limiting case of strong lateral interactions: 
those protofilaments missing at least one neighbor had a very high intrinsic depolymerization rate 
/3qi = 13 s _1 . In this case, model NP produced similar behavior to that seen in experiments only 
with an on rate a full 20-fold higher than our estimate. The maximum depolymerization rate was 
then 2 dimers s _1 or 1 fim min -1 . 

From this analysis we conclude that if kinesin-8 motors depolymerize MTs nonprocessively, 
GMPCPP-stabilized MTs must show some intrinsic instability. Specifically, protofilament inter- 
actions must cause at least one full row of dimers to be destabilized by the removal of a single 
dimer. Based on current experimental results, we cannot unambiguously distinguish this case from 
a processive motor and stable protofilaments. However, we note that the averaged length versus 
time traces are not identical between these two cases. There is less time variation in the MT 
depolymerization rate in model NP than in model P. 

The apparent motor processivity may be due to the inherent instability of MT protofilaments 
lacking neighbors. However, our results indicate that Kip3p likely depolymerizes MTs processively. 



Fluctuations during processive depolymerization 

Our results so far have focused on averaged MT behavior. However, fluctuations of MT length 
about the average are sensitive to kinesin-8 motor depolymerization processivity. To illustrate this 
effect, we chose parameters that give similar average length versus time dynamics by keeping the 
product of the microscopic depolymerization rate fc® n< ^ and the on rate constant k ori approximately 
fixed while varying the maximum depolymerization processivity, a (fig. [3|). When a is increased, 
the dynamics of MT length as a function of time become more rough (fig. [3]A., inset). Motivated 
by the work of Shaevitz et al, who considered the variance of stepping behavior for conventional 
kinesin and showed that the variance increased as the step size increased (j27l ). we quantified the 
fluctuations in MT length as a function of time by determining the standard deviation of MT 
length in our simulations (figure [3j3) . For independent protofilaments, we found that the standard 
deviation increases with a. As expected, the maximum standard deviation scales as a 1 / 2 . For these 
simulations in particular, the maximum standard deviation of MT length ~ 6a 1 / 2 . 

Thus, we propose that experimental measurement of the variance in MT length as a function of 
time can be used to assess depolymerization processivity. One experimental technique to measure 
the fluctuations is to start with MTs of a certain average length, measure the MT length again 
after a fixed time, and determine changes in the width of the MT length distribution. Suppose the 
experimental parameters correspond to those used in fig. [3]A,B. Our idealized simulation started 
with 10-/im long MTs; we then determined the length distribution of MTs after approximately 
1000 seconds of depolymerization, when the variance is largest. The standard deviation of MT 
lengths in the simulations was approximately 320 nm for a = 52, 700 nm for a = 197, and 950 
nm for a = 394, a difference that could be measurable by high-resolution light microscopy. A real 
experiment would begin with a distribution of MT lengths, but would still observe broadening of 
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the MT length distribution. (Note that the precise value of the standard deviation depends on the 
parameters used in the model. Improved measurement of parameters such as the motor binding 
rate would allow more accurate prediction of the fluctuations.) 

In addition, the magnitude of the fluctuations depends on the strength of lateral bonds between 
protofilaments in the MT. We considered simulations without and with strong lateral interactions 
(Po,i = 13 s _1 ), as shown in fig. [3p,D. 

Mean-field model of motor-MT dynamics 

The Monte Carlo simulations accurately represent the full model, but the equations are complex and 
analysis requires running many simulations. We therefore developed simpler mean-field models that 
describe the average occupancy of motors along the MT and the position of the MT end. A mean- 
field model considers averaged values of variables such as motor occupancy and depolymerization 
rate, and therefore neglects stochastic fluctuations. While these models are approximations to the 
full model, they are useful for a simpler, approximate analysis and for understanding the importance 
of fluctuations (which are absent in the mean- field models). 

We analyzed two classes of mean-field model: (i) a description of the motor occupancy profile, 
which characterizes the distribution of motors along the MT, and (it) a description of the MT 
end, which uses the results of the motor occupancy studies to develop a simple model of the MT 
depolymerization dynamics. 

Motor occupancy profile 

To develop and analyze this model, we neglect MT dynamics. (Recall that in the in vitro exper- 
iments, the MT intrinsic dynamics primarily affect the long time behavior, i.e., when the tubulin 
concentration in the flow chamber is sufficiently high to cause significant polymerization.) The 
average fractional occupancy of motors along the MT, p, is described by (|19l ) 

dp dp r /-i \ i /-.\ 
dt = ~ V dx + onC( ~ P '~ oSP ' V) 

On the right hand side, the first term represents the rate of change of motor concentration due to 
biased motion of the motors with velocity v, the second term represents the binding of motors to 
unoccupied sites at rate fc on c, and the third term represents unbinding of motors from occupied 
sites at rate k Q. The bulk motor concentration c is assumed constant. This equation treats motor 
crowding effects in a mean-field approximation: the rate of binding of motors to the MT is assumed 
proportional to (l — p), decreasing in proportion to the average occupancy of a particular site. We 
neglect crowding effects in the transport term ~ , which if considered make the density equation 
nonlinear (fl~9T j. 

The steady-state density distribution away from either of the MT ends is given by the constant 
solution to this equation: 
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Note that if the on rate is sufficiently small, that is, k on c <C k g, then po ~ k nc/k ff = c/K. In 
other words, the average motor occupancy on the MT is approximately the bulk motor concentration 
divided by the equilibrium constant for motor binding to the MT. 

Away from either of the MT ends, the density approaches the constant value pq. If we consider 
a spatially constant occupancy which is not equal to po, the time dependence of eqn. ([I]) has 
exponential solutions. If at time t = motors are introduced to the system, the density far from 
the MT ends will change in time according to 



The characteristic time scale is 



p(t)=p (l-e-^)- (3) 

k oS + k on c' 4 

As shown in figured! this expression agrees well with the simulation results, giving a value of the 
occupancy within 10% of the value from simulations for the region away from the plus end of the 
MT. 

Near the MT ends there is a "boundary layer" where transport effects and boundary conditions 



change the motor density away from po ( la . 1 19 . |20j). Near the minus end (small x) the boundary 



layer has low motor density; this occurs because the number of motors that have moved out of a 
region near the minus end is not balanced by motors moving into that region from the minus end. 
At x = the motor density is exactly zero. In the linear mean-field approximation described by 
equation ([2]), the steady-state density is then 

p(x) = p (l-e-^), (5) 

where A = v/(k g + k on c) = vt is the length scale that characterizes the boundary layer of the 
density distribution near the minus end. In figure (H we show that the steady-state simulation 
results are well represented by this occupancy profile. 

For small on rate (k on c <C k Q g) this length scale is A v/k Q g = £, the motor run length. For 
i < A, the motor occupancy increases linearly with x, with slope po/X = k on c/v. 

Motor occupancy profiles consistent with this result have been observed both in vitro and in 
vivo. Varga et al. observed a linear increase in Kip3p motor occupancy with position in vitro 
(fl5l ). (We note that if the slope of this linear increase could be measured it would provide a direct 
measure of the motor on rate k 0I1 .) Stumpff et al. imaged human kinesin-8 fluorescence along MTs 
in fixed cells and observed a gradient in fluorescence which is qualitatively consistent with the 



model predictions ([13 h similar gradients in motor fluorescence were seen in yeast cells by Varga et 
al. and Gupta et al. (|12l . Il5l). 



Mean-field model of plus-end dynamics 

The dynamics of the microtubule are controlled by the density of motors at the MT end. Here 
we formulate and analyze a mean-field description of the motor occupancy at the MT end and the 
MT length. Because motors promote MT shortening, the density at the MT plus end is constantly 
changing, making the dynamics here typically non-steady state. 
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We define p e (t) to be the average motor occupancy at the last site on the MT end. In this 
model we do not consider protofilament interactions, so we are effectively considering a single- 
protofilament MT. The MT length is represented by L(t) and its rate of change is dL/dt, which 
is negative during depolymerization. The coupled dynamics of the end occupancy and MT end 
change according to: 

(v-^jp(L-e,t)(l-p e )-klf Pe , (6) 

-k^pe. (7) 

The first term on the right side of eqn. ([6]) represents the arrival of motors from the region adjacent 
to the end, where the density is p(L — e, t) and e is a small parameter. Solving for the full time 
dependence of p(L — e, t) is, in general, difficult; one would have to solve for the time varying density 
near a moving boundary. However, we have assumed that a motor falls off while depolymerizing if 
the dimer directly adjacent is occupied (fig. [IP). Thus, if a transient clump of motors developed at 
the MT plus end, it would be quickly removed at rate fc!i nc ^ = 13 s _1 , faster than other processes 
in the model. Thus, we can approximate p{L — e, t) ~ p(L — e) ~ p(L), where p(x) is the motor 
occupancy for a region far from the MT plus end. This density p{x) may vary in time or be a 
steady-state value, but we assume that p{x) is controlled by dynamics away from the MT end. The 
second term on the right side of eqn. ([6|) describes unbinding of the motor at the end. In eqn. ([7]), 
we assume that the rate of MT shortening is proportional to the depolymerization rate and the 
motor density at the end. Note that if p e is constant in time, then according to eqn. (|7|) the MT 
shortens at a constant rate. 

These equations can be combined to write 

^ = (v + k™ d Pe ) p(L)(l - p e ) - k%fp e . (8) 

Equation ([SD can be numerically integrated or studied analytically within certain limits. Below 
we first determine analytically the constant depolymerization rate of very long MTs, where p(L) as 
po as determined from eqn. ([2]). We then compare this predicted constant depolymerization rate 
with simulations on long MTs. 

Then we will consider how long it takes to approach this constant depolymerization rate, and 
find a typical time scale of tens of seconds. This rapid approach to constant depolymerization led 
us to consider quasi-static depolymerization. We will assume that even when the motor density 
away from the MT end varies spatially, the motor occupancy at the MT plus end rapidly tracks 
these changes. 

Constant depolymerization of long MTs 

We will consider first the limit of long MTs with a constant motor density, so that p(L) = po, 
independent of MT length. In this limit, the constant depolymerization velocity of the MT is 
determined by the steady-state value of p e . In this case, equation ([8]) simplifies to one with no L 



dp e 
dt 
dL 
~dt 
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dependence: 

( v + k™ d p e )p (l- Pe )-k™ d p e , (9) 



dp e 

dt 



-k™ d p pl + [(k™ d - v)po - klf]p e + vp . (10) 



This equation has steady state solutions determined by the quadratic equation 

k™ d p p 2 e + [(V - k™ d ) P0 + klf] Pe - V P0 = 0. (11) 

Defining g = (y — k^ d )po + k^^ we can write the solutions as 



-g ± J g 2 + Avk e ^ d pl 
Pe± = • ( 12 ) 

The physically relevant solution (with < p e < 1) is the negative root. When the end occupancy 
is constant, the depolymerization rate is \dL/dt\ = k e ^ ld p e . In fig. [5]we show the predicted steady- 
state occupancy at the MT end, the resulting shortening rate, and comparison with simulations. 
(For details of the simulations, see the section on length-dependent depolymerization below.) 

The effect of fluctuations is to decrease the depolymerization rate, relative to mean-field pre- 
dictions. This is intuitively reasonable, because when a fluctuation leads to a higher than average 
density, the motor at the MT end is then rapidly knocked off (see fig. [1]) , which leads to a decrease 
in the density at the end. Therefore fluctuations which decrease the motor density at the end 
decrease the depolymerization rate more than fluctuations which increase the motor density at the 
end increase the depolymerization rate. Despite this small error, the mean-field theory predicts the 
shape of the curve correctly and determines the depolymerization rate to within 50%. 

Approach to steady state 

We note that these predictions assume that the MTs start with an initial length long enough 
that the steady-state motor occupancy at the end can be reached. This assumption may not 
apply in experiments (see below). We therefore determine the approach of solutions of eqn. (|10p 
to steady state. This equation is a Ricatti equation with constant coefficients, of the form p = 
fp 2 + gp + h, which can be transformed into a linear, 2nd order ODE using the substitution 
u(t) = exp[— j fp(t)dt]. The function u(t) then has two exponential solutions u(t) ~ e r±t with 
inverse time constants 

r ± = l -(g±^/g 2 -Afh). (13) 
The resulting time-dependent solution for the density, given p e (t = 0) = pi, is 



1 fr+(r- + kl nd popi)e r + t -r-(r + + k end pop i )e r - t \ 

pS) = , ot ,a — ; t^a r— — ; 77^a r~rz • ( 14 ) 



fc end pQ I ( r _ + feend pop .) e r+t _ ( r+ + ^.end 



PoPi)e r ~ t 



Typical values of parameters give h < and therefore r + > and r_ < 0, so the dynamics will 
be controlled by the r + terms for long times. The decay times r^} are of order seconds to tens of 
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seconds, with slower times for lower concentrations. See fig. [6] for the values of the time constants 
and a typical time trace determined by eqn. (|14p . 

Note that this equation assumes that the motor density outside the end of the MT has reached 
the steady-state value pQ. As shown in fig. [U it can take several minutes for the density (away from 
the MT end) to equilibrate. Therefore the dynamics of the MT end are limited by the dynamics 
of motors away from the end more than by the processes at the end. 

Quasi-static depolymerization of MTs 

Outside the regime of constant depolymerization, we can still make a simple approximation to the 
depolymerization rate by assuming that the plus end of the MT has a motor density determined by 
the instantaneous solution of the steady-state equation, but with a varying density away from the 
end. In other words, we assume that the time for motor density at the MT plus end to reach steady 
state is short compared to other time scales in the problem. This is a reasonable approximation, 
given that the dynamics at the end of the MT reach steady state in seconds to tens of seconds 
(fig. [6]) while MT shortening typically takes minutes. 
In the quasi-static approximation, we solve 



dL 



= [v-^jp(L)(l- Pe )-klfp e , (15) 

§ = -*- n v < 16 > 

The solutions here are similar to those of eqn. (j 12 [) . but with the varying density p(L) apparent in 
the solution: 

k^p{L)pl + [(v - k™ d )p(L) + klf]p e - vp{L) = 0. (17) 
Defining g(L) = (v — k^^)p(L) + k^^~, the physically relevant solution is 



-g(L) - x /g(L)2 + 4v k™ d p(Ly 
2k™ d p{L) 



This quasi-static density at the end of the MT then determines the depolymerization rate via 
eqn. (|16p . We will study the accuracy of this quasi-static approximation below in the section on 
the phase diagram of length-dependent depolymerization. 



Length-dependent depolymerization 

A key feature of the Varga et al. results is the observation of length-dependent MT depolymer- 
ization (fl5l). where the depolymerization rate decreases as the MT length decreases. Since the 
depolymerization rate cannot increase indefinitely with MT length, a long MT (L 3> €) will not 
show length-dependent depolymerization until it becomes sufficiently short. The "crossover" length 
d is the length at which length-dependent depolymerization begins; for MT lengths L > d the de- 
polymerization rate is constant while for L < d the depolymerization rate decreases as L decreases. 
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Here we determine the crossover length d as a function of motor parameters and experimental 
conditions. 

Motors accumulate at the MT plus end for several reasons. First, motors have different residence 
times at the MT end than away from the end. Away from the ends of the MT, a motor moves away 
from a dimer at rate v. The off rate for a motor from the last dimer along the MT is k®^. This 

will increase the motor occupancy by a factor of approximately v/k e ^ ~ 200 (although note that 
the average motor occupancy per site cannot be > 1; see fig. [5]). Second, depolymerization moves 
the MT end closer to motors away from the end; in the frame of the MT end motors approach at a 
rate v — dL/dt, where dL/dt is the rate of change of MT length. Third, motors could accumulate 
at the MT end due to direct binding of motors to the end. The latter effect is neglected here. 

To characterize length-dependent depolymerization and to understand its role in vivo, we es- 
timated the crossover length d analytically and compare it to simulations. We consider length- 
dependent depolymerization in two limits: first we consider the case most independent of initial 
conditions, which corresponds to starting with very long MTs or pre-equilibrated motor occupancy 
along the MT. Second, we consider the case where the initial motor occupancy on the MT is 
zero; this case is more relevant to an experiment started with the addition of motors to previ- 
ously unoccupied MTs. See fig. [7] for a comparison of depolymerization of pre-equilibrated MTs to 
depolymerization with an initial transient. 

Length-dependent depolymerization independent of initial conditions 

Here we characterize length-dependent depolymerization in a regime which is, as much as possible, 
independent of the initial conditions or the starting time of an experiment. Such a limit would be 
reached in experiments (or in cellular conditions) if one either (£) started the experiment with very 
long MTs, or («) equilibrated the motor density on the MTs before the start of depolymerization. 

To approach this limit in the simulations, we implemented both long initial MTs and motor 
pre-equilibration. We started with MTs 4000 dimers long (32 fj,m) and ran the simulations for 1000 
seconds with motor binding, unbinding, and motion allowed but depolymerization turned off. This 
starting point was then used for simulations of MT shortening. Averaged simulation results are 
shown in fig.[8j By examining the depolymerization rate versus MT length (fig. ©3), we can see that 
initially the depolymerization rate is high, but it drops quickly as the clump of motors at the MT 
plus end is removed during depolymerization. Then the depolymerization rate reaches a constant 
value that holds over a range of MT lengths. (This constant rate is compared to the predictions 
of the mean-field model in fig. [5j) Once the MT is sufficiently short, the depolymerization rate 
drops below the steady-state value. Eventually the depolymerization rate drops to zero for a MT 
of length zero. 

We also determined the fluctuations in MT length about the average: fig.EP shows the standard 
deviation in MT length as a function of MT length during shortening. Lower motor concentration 
leads to a significantly larger standard deviation in MT length, indicating that larger fluctuations 
about the average occur when fewer motors are present. 

The length where the depolymerization rate drops below the steady-state value is the "crossover 
length" d where length-dependent depolymerization sets in (fig. [Hp). We defined the crossover 
length d as the MT length where the depolymerization rate first decreased by 20% from the steady- 
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state value. The results from simulations are shown in fig. [Hp. We calculated d in the mean-field 
model by rewriting eq. ([HJ using L rather than t as the independent variable: 

*~ = (v + k™ d p e )p(L)(l-p e )-klf Pe , (19) 



dp e 

dL fcend 



Pe 



v + k™ d p e ) p(L)(l - p e ) - k e Q fp e ] . (20) 

We then numerically integrated eqn. (|20p with Lq = 2 x 10 5 dimers and determined the crossover 
length d where the depolymerization rate decreases by 20% below the steady-state rate. These 
results are the solid curve in fig. [Hp. The predictions of the mean-field model agree well with the 
simulations. 

For these results, we find that the crossover length d is 10 pm or longer for bulk motor con- 
centrations of 10 nM or lower. However, for high bulk motor concentration of 50 nM or more, the 
crossover length decreases to 4 pm or less. This suggests that length-dependent depolymerization 
is only prominent for sufficiently low motor concentration. This difference may partially explain 
why Varga et al. observed length-dependent depolymerization but Gupta et al. did not. The Varga 



experiments used lower bulk motor concentrations (|12|, [15|). 

To further understand the behavior of length-dependent depolymerization, we examined how 
the crossover length d changes when two key parameters are varied. We changed the bulk motor 
concentration and the motor off rate at the MT end and determined d. We studied this crossover 
using the mean-field model in two limits: first, we numerically integrated eqn. (|20p to determine the 
crossover length in the full mean-field model. Second, we determined the crossover length in the 
quasi-static model, which assumes that the motor occupancy at the MT end is instantaneously de- 
termined by the solution to the steady-state eqn. (118p . This solution is valid if the motor occupancy 
at the end changes quickly compared to other time scales in the problem. The results are shown 
in fig. GO We note that the quasi-static approximation gives results similar to the full mean-field 
model except at the lowest motor off rates. This is intuitively reasonable, since slower motor off 
rate at the MT end increases the time scale for motor dynamics at the MT end, decreasing the 
validity of the quasi-static approximation. 



Length-dependent depolymerization controlled by initial conditions 

The experiments of Varga et al. do not begin with very long MTs such as those discussed in 
the previous section. As a result, we expect the dynamics to show a lag relative to the steady 
depolymerization observed when starting with long MTs (see fig. [7]). To address conditions relevant 
to experiments, in this section we characterize the effects observed in this transient regime. 

For the mean-field calculations presented here, we numerically integrated eqn. using a time- 
dependent motor density away from the end: 

^ = {v-^)p(L,t)(l-Pe)-k e fp e , (21) 
dL , PTir i . . 
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The initial conditions are L(t = 0) = Lq and p e (t = 0) = 0. The time dependence of the density 
away from the end, p(L, t) is described by 

p(x,t)= Po (l-e- x / x )(l-e- t / T ), (23) 

with r = l/(k Q + k on c) and A = vt (as discussed above; see eqns. (|3|), ([!]), and fig-H]). This form of 
the density is an approximation that assumes the density approaches its steady-state distribution 
with dynamics controlled by the slowest time scale in eqn. ([2]). Note that we assume that at time 
t = motors are introduced to the system, as in the Varga et al. experiments. 

In fig. [10] we illustrate the solutions to these equations for an assumed bulk motor concen- 
tration of 5 nM. The dynamics vary with the initial length of the MT assumed; the maximum 
depolymerization rate varies by almost a factor of 10 with a comparable change in Lq. 

We compared these results to the experimental data of Varga et al. by considering varying 
initial MT lengths and bulk motor concentrations. For each curve, we determined the maximum 
depolymerization rate (the peak of the curve in fig. [TUB). The results are shown in fig. [TTj where 
we show how the maximum depolymerization rate varies with Lq and bulk motor concentration. 
The slope of depolymerization rate versus initial MT length is shown as well, with the results of 
Varga et al. shown for comparison. The results are reasonably similar — the model predictions are 
within a factor of 2 of the experimental results. 

The model predictions show a linear dependence of the slope of depolymerization rate on initial 
MT length, while the Varga et al. data suggest a non-linear dependence of the slope on bulk motor 
concentration (fig. [UJ). This nonlinearity could result from motor cooperativity, which would lead 
to a nonlinear dependence of the depolymerization rate on the motor density at the MT end. This 
would be an interesting direction for future experiments to explore. 



Conclusion 

We have developed a theory of MT depolymerization by kinesin 8 and demonstrated agreement 
between our theory and currently available experiments. The model incorporates biased motor 
motion toward MT plus ends, motor-catalyzed depolymerization of plus ends, motor binding and 
unbinding, and motor crowding effects. Our theory quantitatively reproduces the experiments of 
Varga et al. using purified Kip3p, the budding yeast kinesin 8. 

In developing a theory of kinesin-8 motors, we addressed a limit which has not been considered 



in previous theoretical work (|la . 1 19 . [2C], |22j). Experiments have revealed that clumps of kinesin 



8 form on the MT plus end and change significantly during MT depolymerization (|12l . Il5l ). This 
observation suggests that MT shortening and motor density changes occur on similar time scales. 
Therefore a steady-state mathematical analysis is likely to be poor approximation to experiments. 
We have therefore developed and analyzed a time-dependent equation for the MT end that couples 
MT depolymerization, motor arrival and dissociation, and motor crowding effects. The results of 
the equation for MT end dynamics agree well with full Monte Carlo simulations of the model. 

Despite the good agreement between our theoretical results and experiments, there are several 
uncertainties in our model. We have estimated the motor binding rate constant, which has not been 
precisely quantified in experiments. The most significant effect of this uncertainty is that it leads 
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to uncertainty in the depolymerization processivity, as illustrated in fig. [3l In addition, changing 
the on rate constant would alter the dependence of the crossover length d on motor concentration, 
changing the predictions of fig. [U We also assumed that when a motor catalyzes dimer removal, it 
falls off the MT if another motor is bound directly behind it. This assumption is supported by in 
vivo experiments with fluorescently-labeled Kip3p, which show that a clump of motors accumulates 
at the MT plus end when the MT is growing, but the clump shrinks during MT shortening (|12l . Il5l ) . 
Finally, we did not consider possible effects of motor cooperativity in depolymerization; cooper- 
ativity is a possible explanation for the apparent nonlinearity of the experimental data shown in 
fig- [HI Our best-fit parameters did lead to quantitative agreement between the theory and the 
experiments of Varga et al. Errors in our determination of the parameters would lead to changes 
in the quantitative predictions of the theory, but we verified that the qualitative predictions are 
insensitive to the exact parameter values chosen. 

We observed length-dependent depolymerization of MTs in the model, as seen experimentally 
by Varga et al. (fl5l ) and theoretically by Govindan et al. (fl8l ). However, this phenomenon occurs 
only below a "crossover" length d; MTs much longer than d will depolymerize at a constant rate. 
The crossover length is controlled by key parameters of the model, particularly the bulk motor 
concentration, the motor translocation processivity, and the depolymerization processivity. In 
particular, using the best-fit experimental parameters we predict that the crossover length will 
decrease from 16 //m to 4 /im if the bulk motor concentration is increased from 5 to 50 nM. This 
strong concentration dependence may explain the observation of length-dependent depolymerization 
in vitro by Varga et al. (fl5l ) and not by Gupta et al. (12j). (Note that the experiments differed in 
other ways, such as the observed motor velocity.) 

Length-dependent depolymerization is not specific to motors with biased motility (fl8| ). For 
example, length-dependent depolymerization could occur for the kinesin-13 MCAK, which diffuses 
on MTs. The crossover to length-dependent depolymerization is sensitive to motor processivity. As 
a result, we expect length-dependent depolymerization will be much more important in the kinesin- 
8 family of proteins because they have a much higher processivity than do kinesin 13s (~ 12 ^m 
for Kip3p vs ~ 1 /xm for MCAK) (0). 

The fact that the crossover length is so parameter dependent highlights the need for care in 
interpreting experimental results. In cells, the bulk motor concentration and other parameters 
may be different from the in vitro values. If the key features of our model do apply in the more 
complicated cellular environment, we can make some speculative predictions. In particular, an 
overexpression experiment, which increases motor concentration, would lead to both shorter MTs 
and to a decrease in the crossover length d, and therefore a decrease in the range of lengths over 
which there is length-dependent depolymerization. 

Our results are most consistent with processive depolymerization by kinesin-8 motors, suggest- 
ing that multiple tubulin dimers can be removed by a single motor. In particular, the experimental 
data can be fit to a non-processive model only if (i) the motor on rate is a factor of 20 higher than 
the one estimated from experiments, and (ii) GMPCPP MTs are intrinsically unstable, so removal 
of one dimer from one row of tubulin subunits leads to the removal of all 13 dimers around the 
MT. We note that processive depolymerization has been observed for the kinesin-13 motor MCAK 
(j2ll ). so kinesin 13s and 8s may share some features in their depolymerization activity. 

We have shown that processive depolymerization tends to increase the fluctuations of MT length 
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about its average (fig. [3]). Therefore the roughness of the MT length versus time during depolymer- 
ization could be used to quantify the depolymerization processivity. While we emphasize that our 
theory may not be directly relevant in cells, our results on MT length fluctuations do have inter- 
esting implications for some recent experiments. Although kinesin 8s serve, on average, to decrease 
MT length, the large fluctuations that result from processive depolymerization would introduce a 
significant variance about this average behavior. In addition, higher protein concentrations should 
decrease these fluctuations (fig. [DC). If this were true in cells, it could explain two recent and puz- 
zling observations on kinesin 8 in vivo. Stumpff et al. found that kinesin 8 over-expression decreased 
the amplitude of metaphase chromosome oscillations, while reduction in kinesin-8 concentration by 
RNAi increased it (fl3l). Unsworth et al. found that deleting the kinesin-8 genes (completely elimi- 
nating the protein) from fission yeast decreased MT length fluctuations ()14l ). as would be expected 
if physiological levels of kinesin 8 contribute significantly to length fluctuations. Although our 
simple mathematical model does not describe the full complexity of the mitotic spindle, it does 
highlight the potential role of motor-induced fluctuations in spindle behavior. 
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Table 1. 



Quantity 


Symbol 


Typical value(s) 


T~» • 1 ' 1 TV Tm 

Position along Ml 


X 


- 2500 dimers 


MT length 


L 


- 2500 dimers 


Motor velocity 


V 


7.5 dimers s _1 


Motor on rate constant 




3x10 site 1 nM 1 s 1 


Motor off rate 


k oS 


5 x 10" 3 s" 1 


Equilibrium constant for motor 
binding to MT 


K = /c n «7 /c n 

Ull ' v - ,j -- L 


1.67 fiM 


Motor run length 


e = v/k oS 


1500 dimers 


Motor off rate at MT end 


i,end 


3.3 x 10~ 2 s" 1 


Rate of motor-catalyzed MT de- 
polymerization 


£,end 


13 Himprs s — 1 


Motor depolymerization proces- 
sivity (upper bound) 


a = k end /k er ^ d 
I off 


1 — 400 dimers 


Bulk motor concentration 


c 


1 - 200 nM 


Motor occupancy per tubulin 
dimer 


p(x) 


0-1 


Steady-state motor occupancy 
away from MT ends 


PO = k on c/ (fc on c + k oS ) 


6 x 10~ 4 - 6 x 10~ 2 


Motor occupancy at MT plus end 


Pe 


- 1 


Motor occupancy boundary 
length 


X = v/ (k on c + k oS ) 


1400 - 1500 dimers 


Time scale of approach to steady- 
state motor occupancy 


t = 1/ (Vc + k oS ) 


185 - 200 seconds 


Length of crossover to length- 
dependent depolymerization 


d 


- 2400 dimers 


Bulk tubulin concentration 




10 - 100 nM 


MT polymerization rate constant 


a 


5.4 dimers ^M _1 s _1 


MT depolymerization rate 


P 


0.1 dimers s _1 
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Figure Legends 

Figure [TJ 

Model of kinesin-8 motor protein's interaction with a MT protofilament showing the key rates. 
Rates in black affect only the motor, while those in red affect the MT plus end. The plus end 
of the MT protofilament is indicated by a thick vertical (red) line, and the dimers are indicated 
by (blue) boxes. Note that depolymerization (at rate k^ 1 ^) affects both the motor and the MT 
plus end. (A) A motor binds to a dimer of the MT with on rate fconC and unbinds with rate k Q g. 
The motor steps forward at rate v; backward motion is not considered due to the biased motion of 
kinesin 8s. MT dynamics are represented by allowing dimers to add to a MT end at rate act (where 
ct is the bulk concentration of tubulin dimers) and dissociate at rate (3. (B-D) Depolymerization 
models. (B) If the motor depolymerizes processively, it removes a MT dimer as it steps backward 
(with rate k^ 1 ^), thereby shortening the MT. (C) If the dimer behind the MT end is occupied, 
the motor falls off the MT in either model. (D) In the nonprocessive depolymerization model, the 
motor removes a single tubulin dimer and falls off the MT. (E) Lateral interactions help stabilize 
MTs. We incorporated this into our model by allowing the depolymerization rate to depend on the 
number of neighboring protofilaments (i.e., 0, 1, or 2). In this case, the rate at which a terminal 
tubulin dimer unbinds from a protofilament is given by (3q if the dimer has no lateral neighbors, (3\ 
if the dimer has 1 lateral neighbor, and fa if the dimer has 2 lateral neighbors. 



Figure El 

The slow polymerization of GMPCPP MTs in the presence of free tubulin accounts for the long- 
time tails observed in Varga et al. (|15l ). In that work, Kip3p was unable to fully depolymerize 
the MTs over the course of a single experiment. In our simulations, adding previously measured 
MT polymerization and depolymerization rates for GMPCPP stabilized MTs (0) reproduced the 
observed behavior. The solid traces were made assuming that the only MT dynamics were those 
caused by the motors, while the dashed traces were made including intrinsic MT polymerization 
and depolymerization. 



Figure G3 

Effects of depolymerization processivity on MT length fluctuations. We chose model parameters 
which led to similar overall behavior of MT length versus time (A,C) but had different motor- 
induced depolymerization rates and motor on rates. All simulations of model P used the experi- 
mentally derived value k^^ = 3.3 x 10~ 2 s _1 . The maximum processivity (the maximum numbers 

of dimers removed per motor) is a = k^^/k^^. (A) We first considered completely stable MTs, 

those with /?o,i,2 = 0, and varied k on and fc® to obtain curves which all have a similar average 
shape. The trace with a = 394 uses c = 1 nM and the best-fit parameters found when comparing to 
experiments: k on = 3 x 1CP 6 site -1 nM" 1 s -1 and k 611 ^ = 13 dimers s —1 . The curve with a = 197 
has the on rate constant doubled to 6 x 10~ 6 site -1 nM -1 s -1 and the maximum depolymerization 
rate halved to 6.5 dimers s . The curve with a = 52 has the on rate constant increased by a factor 
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of 8 to k n = 24 x 10 site nM s and k^ 1 = 1.7 dimers s . Curves in the main panel 
show averages of 500 simulations, each of a MT with 13 independent protofilaments. Curves in 
the inset panel show results of individual simulations; the roughness of the MT length versus time 
behavior decreases as a is decreased. (B) The standard deviation of MT length as a function of 
time for the simulations shown in (A). Simulations with higher a show larger standard deviation. 
(C) Allowing protofilaments without two neighbors to spontaneously depolymerize (P, A; n, high 
0) gives a similar average curve to a simulation with fully stable MTs, but fc n increased by a 
factor of 13. A similar average curve can also be obtained in the non-processive case, but only if 
protofilaments without two neighbors spontaneously depolymerize and k on is increased by a factor 
of 300 (NP, 300 x A; on , high (3). (D) The standard deviation of MT length for the simulations shown 
in (C). 

Figure [4l 

(A) Examples of motor occupancy profiles. For this figure, we assumed a much higher on rate 
than found in experiments to make occupancy changes visible. The reference parameter set has 
konC = 0.002 dimer -1 s -1 (which would correspond to a bulk motor concentration of 667 nM at the 
typical on rate constant of k on = 3 x 10~ 6 dimer -1 s -1 nM -1 ), k Q = 0.005 s -1 , and fc^jj = 0.02 
s -1 (black solid line). Therefore I = 1071 dimers (black vertical line). The curve with decreased I 
(blue dotted line) has feonc and k g both doubled to halve t to 536 dimers (blue vertical dashed- 
dotted line) while keeping Krj unchanged. For these parameters, the mean-field expression for 
the occupancy from eqn. ([5]) is the blue dashed/dotted line. The curve with decreased Kjj has 
k nc halved (red dashed line). (B) Average motor density as a function of time for the lower £ 
parameter set at three positions along the MT: £/2 (black solid line), 3^/2 (red solid line) with the 
mean- field analytic expression of eqn. ([3]) superimposed (black dotted line), and the MT plus end 
(blue dashed-dotted line). All plots are averages of 500 simulated MTs. 

Figure [5l 

Steady-state motor occupancy and depolymerization rate of long MTs. Left axis shows steady- 
state motor occupancy at the MT plus end (solid blue line, black circles) or away from the end (red 
dashed-dotted line) as a function of bulk motor concentration. The mean-field model (solid blue 
line) is the prediction of equation (fT2"j) for the steady-state occupancy of the MT plus end. Right 
axis shows the resulting steady depolymerization rate in the mean-field model (solid blue line) and 
simulations (black circles). The simulation results were determined from the simulations shown in 

fig. El 

Figure O 

Approach to steady-state end occupancy in the mean-field model. (A) End occupancy, beginning 
from p e (t = 0) = 0, for representative values of the bulk motor concentration. (B) Values of the 
time constants r^ 1 as a function of bulk motor concentration. 
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Figure 

Steady-state versus transient depolymerization. The transient condition (when motors are added at 
t = 0, red dashed line) gives dynamics with an initial lag when compared to the condition with pre- 
equilibrated motors (blue solid line). Results shown are from simulations (average of 500 runs) with 
bulk motor concentration of 5.5 nM. (A) MT length as a function of time. (B) Depolymerization 
rate as a function of MT length. 

Figure [8l 

Depolymerization of long MTs. Each curve is the average of 500 independent simulations. Each 
simulation was started from a pre-equilibrated MT: the simulation was run for 1000 sec with no 
filament depolymerization, to allow the motor density on the MT to reach steady state. (A) Length 
versus time. (B) Depolymerization rate versus MT length. Black squares indicate the crossover to 
length-dependent depolymerization. (C) Standard deviation of MT length versus MT length. (D) 
Length of crossover to length-dependent depolymerization in the simulations (red circles) and the 
mean- field model (blue solid line). 

Figure [|0 

Length of crossover to length-dependent depolymerization, as a function of the bulk motor concen- 
tration and the motor off rate at the MT end. The horizontal dashed line is the motor off rate at 
the MT end estimated from the experiments of Varga et al. (A) Mean-field model. (B) Quasi-static 
approximation, where the motor occupancy at the MT end is assumed to change quickly compared 
to other dynamics in the problem. 

Figure I10L 

Dynamics under conditions of transient shortening depend strongly on the initial MT length. Black 
solid lines: simulations (average of 500 runs). Blue dashed lines: mean-field model. At t = 0, motors 
are introduced to the system. The bulk motor concentration was 5 nM. (A) MT length as a function 
of time. (B) Depolymerization rate as a function of MT length. 

Figure II 1L 

Variation of the maximum depolymerization rate with initial MT length, showing predictions of 
(A) the simulations and (B) the mean-field model under conditions of transient shortening. (C) To 
compare to the results of Varga et al., we fit a straight line to the first 1000 dimers of each curve 
in (A-BV and determined the slope as a function of bulk motor concentration. The data of Varga 
et al. (|15l ) are shown for comparison. 
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